Optimizing the classical heat engine. 
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A pair of systems at different temperatures is a classic environment for a heat engine, which 
produces work during the relaxation to a common equilibrium. It is generally believed that a direct 
interaction between the two systems will always decrease the amount of the obtainable work, due 
to inevitable dissipation. Here a situation is reported where, in some time window, work can be 
gained due the direct coupling, while dissipation is relevant only for much larger times. Thus, the 
amount of extracted work increases, at the cost of a change of the final state. 
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Introduction. — The classical problem of thermody- 
namics is the determination of the maximal amount of 
work that can be extracted from a non-equilibrium sys- 
tem, during the relaxation to the equilibrium state [Q-Q ■ 
In recent years the interest in this time-honored prob- 
lem was renewed (see [Q for review). In many cases the 
limits proposed by the founders of the classic thermo- 
dynamics have appeared as too idealistic, and attention 
was focused on the study of dissipative effects, which re- 
strict the abilities of realistic heat engines Q . One of the 
most popular examples of that kind is a direct interac- 
tion between the thermal baths that drive the standard 
heat engine || . Evidently, this is a way to dissipate en- 
ergy. The common opinion, expressed in textbooks (ij-jjj 
is that dissipation is the main, if not the only effect of 
any direct interaction. 

The purpose of the present paper is to show that a 
direct interaction between the baths may have energy 
transfer, rather than energy dissipation, as its main phys- 
ical effect. The reason is that the relevant timescale of 
the transfer arising from the direct coupling can be widely 
separated from the dissipative timescale. This leads to 
an optimization of the classical heat engine, in a certain 
time window. As a result, the final equilibrium state of 
the heat baths is changed. The mechanism arises from 
the recently proposed steady adiabatic state Q (see also 
H0 in the related context). 

First, we shall reproduce the classical discussion of the 
maximal amount of work that can be extracted from a 
non-equilibrium system. The underlying ideas are, of 
course, well-known, and we shall need them for a careful 
interpretation of our main result. 

The classical analysis. — Consider two subsystems 
which have different temperatures Ti and T 2 . Depending 
on the concrete relaxation process, the whole thermally 
isolated system will produce different amounts of work. 
This model is well-known and described in textbooks on 
thermodynamics Jl|-§|. We shall assume that the total 
volume of the system is unchanged by the relaxation (al- 
though it can vary during the process), since we wish to 
take into account only the work that can be done due to 
the non-equilibrium initial state, and not due to the gen- 



eral expansion. Denoting by [/,-, Si the initial energy and 
entropy of the system, we get the following expression for 
the work W performed by the system 



W = Ui - U(S), 



(1) 



where U (S) is energy of the final equilibrium state as 
a function of its entropy. Because the temperature is 
positive, U is a monotoneous function of S (dU/dS\v = 
T > in equilibrium). Therefore, W is maximal when S 
is as small as possible. Since the whole system is ther- 
mally isolated, the second law demands S > Si. The 
maximal amount of work is attained for S = Si, i.e., for 
a reversible process towards equilibrium. It is believed 
that for obtaining the maximal amount of work any di- 
rect interaction between subsystems should be removed, 
because it would induce irreversible relaxation. Thus, a 
third body ("engine") should operate between our sub- 
systems to perform the work. At the end of the relaxation 
process the working body must return to its initial state. 
The expression of the total amount W of the extracted 
work can be given as 



W = Ux{Tx) + U 2 {T 2 ) - E7i(T ) - U 2 (T ), 



(2) 



where the final temperature To is determined by the re- 
versibility condition 



Si(Ti) + S 2 (T 2 ) = 5i(T ) + S 2 (T ), 



(3) 



and Uk, Sk (k = 1,2) are energy and entropy of the cor- 
responding subsystems. The reversible relaxation pro- 
cess consists of an infinite amount of elementary cycles 
of the third body. The famous Carnot process is the 
most representative example for such a cycle. The local 
efficiency rj of the optimal cycle, defined as the maximal 
extracted work during the cycle divided by the input en- 
ergy, depends only on the conditions of reversibility and 
conservation of energy: r) = 1 - (Ti/T 2 ), if T x < T 2 |]-|. 
This efficiency is universal and system-independent, em- 
phasizing the power of thermodynamics. In contrast, the 
maximal amount of the extracted work (^|) is not univer- 
sal, and can vary from one system to another. 
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The steady adiabatic state. — Steady adiabatic systems 
have two distinctive properties ||: 1) different subsys- 
tems are not in the mutual equilibrium, but possess dif- 
ferent temperatures, 2) the local characteristic relaxation 
times of these subsystems are very well separated. Let 
us consider a statistical system which has two subsys- 
tems with coordinates x\ and X2- {x\ and X2 can thus 
also code also a set of variables; we will not denote that 
explicitly.) 

The corresponding relaxation times are denoted by Fi 
and r 2 . The condition of well separated timescales is en- 
sured by 7 = Ti/r2 <C 1. The stationary distribution of 
this system can be found from the following heuristic ar- 
guments (for a more rigorous presentation see Q). Due 
to the large difference between relaxation times the x\- 
subsystem comes to an equilibrium while the a^-variable 
is almost unchanged. The corresponding Gibbs station- 
ary distribution reads 



P{xi\x 2 ) 



1 



Z{X 2 



■ exp[-/3 1 H(x 1 ,X2)}, 



(4) 



where 



H(x 1 ,x 2 ) = H 1 (x 1 ) + H 2 (x2) + gH lnt (x 1 ,x 2 ) (5) 

is the total system's Hamiltonian, T\ = 1/ (3\ is the tem- 
perature of the xi-subsystem, and Z{x 2 ) is the partition 
sum at fixed X2- For obtaining the coarse-grained dis- 
tribution P(x2) notice that after integrating out of the 
vast variable x\, the effective Hamiltonian for the slow 
variable X2 is just — Ti In Z(x2), the free energy of the 
Xi-subsystem at fixed x 2 . The stationary distribution of 
X2 then reads 



P(x 2 ) 



Z T ^(x 2 ) 



Z = J dx 2 Z Tl / T2 (x 2 ). 



(6) 



The complete stationary distribution can now be written 
as P(xi,X2) — P(x2)P(xi\x2)- The mean energy and 
entropy of the system are given by the following general 
definitions 



U=f dxi_dx 2 H(xi,x 2 )P(xi,x 2 ), 
S = —J dxidx 2 P(xi,X2)\nP(xi,X2) 



(7) 
(8) 



S can be decomposed as the sum of the entropies of the 
slow and fast subsystems, S = Si + S2, with 



Si = J dx 2 P{x 2 )[- J dx 1 P{x 1 \x2)\nP{x 1 \x 2 )], 



s 2 = - 



J dx 2 P{x 2 )\nP(x 2 ). 



(9) 



Si is the entropy of the fast variable x\ at X2, averaged 
X2 , while S 2 is the entropy of this slow variable itself. As 
was shown in J8| the considered system admits a ther- 
modynamical description. In particular, defining the free 
energy as F = — T 2 InZ, we get: 



F = U — TxSi - T2S2 



(10) 



where the entropies can also be obtained by the standard 
relations: 



Si 



dF 

Wi 



s 2 = - 



dF 
df 2 



(11) 



The energy and entropies are constant in the steady state, 
but a direct coupling induces a steady entropy produc- 
tion at rate S and an energy dissipation at rate II. These 
quantities were analyzed in Ref . [^) , and the obtained for- 
mulas read II = T 2 S + C(7 2 ), with 



S = J 9 - TiT 2 T2 

dH int (xi,x 2 ) _ 
8x2 



dx2 dxi P(xi,X2) x 



dyP(y\x 2 



dH int (y,x 2 ) 



-1 2 



Ox? 



0( 7 2 ) 
(12) 



Although these results were obtained for the strictly 
steady state, they can be applied also in the time- 
dependent case, if the characteristic time of this quasi- 
stationary process is much larger than the largest re- 
laxation time, r 2 . For instance, to obtain the time- 
dependent distribution function for the case of slowly 
(adiabatically) changing temperatures or other parame- 
ters, one just inserts these time-dependent values directly 
in Eqs. ([|, ^|). In this context, the change of free energy 
( |l0| ) can be shown to be the adiabatic work W a d done 
on the system, when varying a parameter a (for exam- 
ple, the width of the potential, or a coupling constant) 
from the initial value to final value a/ at constant 
temperatures: 



r a f r 

Wad = I da dx 1 dx2P(x 1 ,x 2 ,a) 
= F(T 1 ,T 2 ;af)-F{T l ,T 2 ;ai) 



dH(x l7 x 2 ,a) 



da 



(13) 



This is fully analogous to the property of the usual 
(single-temperature) free energy. 

The dissipative effects (given by Eq. (|lj)) are small for 
small 7, g. Let us neglect them for the moment; later we 
shall show that this is allowed in a certain timewindow. 

The maximal amount of work extracted from the steady 
adiabatic state. — Certainly, we can apply the above- 
mentioned general analysis, concerning the maximal 
work, to our adiabatic system. In fact, this analysis 
does not use any concrete property of the initial non- 
equilibrium state, but we should take into account that 
our subsystems interact directly, and not only through 
the third body. As compared to the case without direct 
coupling, the system will now relax to a different equilib- 
rium state, and this is the reason why one can get more 
work done by it. The total amount of the gained work 
can be again written as 
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W(g) = U(T 1 ,T 2 )-U(T g ,T g ), 



(14) 



and the temperature T g of the final equilibrium state is 
defined from the condition of reversibility 



S{T X ,T 2 ) = S{T g ,T g ), 



(15) 



This condition involves the total entropy of the inter- 
acting subsystems, and now assumes that there are no 
additional sources of dissipation besides (fl2|). 

We shall investigate Eqs. ([l4[ |l5|) to first order in the 
small parameter g. Hereafter quantities of the order g° 
and g 1 will be indicated by the index and 1, respec- 
tively. It is evident from (fy that to order g 1 it holds 



F = F (0) + gV 



(16) 



where 



V (Ti,T 2 ) = J dxidx 2 P(xi,X2)H in t(xi,x 2 ). (17) 

Using Eqs. (00), one gets 
Si = Si 0) + gS^ = S{ 0) -g&r.Vo, 



S 2 = sf + gS { 2 1] = Si 0) - gd T2 V , 



(18) 
(19) 



S = + gSW = S<°> - g{d Tl + dr 2 }V , (20) 
U = + glfi 1 ) = + g{l - T t d Tl - T 2 d T2 }V (21) 

Let us now obtain from Eq. (|l^) an expression for the 
final temperature T g to order g, given the value of T , 
the final temperature for 5 = 0, 

T ' = ^{ 1 + s (h + ch )■ (22) 

Here C k = T k d S { ° ] / dT k \ v = dU { ° ] /dT k \ v , with k = 
1,2, are the heat capacities of the subsystems when they 
are uncoupled. Starting from Eq. (|lj) and using Eqs. 
(|l8|)-(pT|), we finally derive the excess work at order g 



(23) 



w{1) = Um W(g)-W(0) 
9^0 g 
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= VoiTuTi) - V (T ,T ) + ^(T - T fc )d Tfc V (Ti, T 2 ). 

k=l 

This is the first main result of our work. It remains to 
be shown in a specific example that this quantity can be 
positive. Let us first point out that further simplifica- 
tions occur when T\ is close to T 2 . To first order in the 
parameter T\ — T 2 eq. (3) gives us 



rp C\ +C 2 

Tq Ci t 1 + c 2 t 2 TiT ^ 

while Eq. (|23|) can be approximated by 



(24) 



2 

E 

k=l 



(To - T k )(d Tk V - {d Tk V }\ Tl =T 2 =T ). (25) 



To illustrate the general results let us present a con- 
crete model, where the direct interaction increases the 
total amount of work: W(g) > W^(0). One of the most 
popular models of the thermal bath is a set of harmonic 
oscillators which is frequently used to derive kinetic 
equations or to gain fundamental insight. Following this 
well-established tradition, we shall model our first (sec- 
ond) thermal bath by Ni(N 2 ) oscillators at temperature 
T\ (T 2 ), and weakly- anharmonic interaction: 



. Ni . N 2 N 



2^ L ' 1 2 

i=l i=l 



i x 2,iJ 



(26) 



where g > 0, and Nx,N 2 > N. It is straightforward to 
show that to order g one has the partial partition sum 

/ , n 2 n \ 

z(x 2 ) = if l/2 exp --hJ2 x h ~aJ2 x h ( 2? ) 

\ i=l i=l / 



and the full partition sum 

Z = T ? NlTl/T2 t} N2 (1 + 2gT x )-^ N 



(28) 



The latter result yields the free energy F = — T 2 In Z 

F = —NiTi In Ti - ~N 2 T 2 In T 2 + gNT x T 2 (29) 

in agreement with the fact that Vb = NT{F 2 . According 
to previous rules we derive 

5 1 = ^N 1 (lnT 1 + l)-gNT 2 (30) 

5 2 = ^N 2 {\nT 2 + 1) - gNTi (31) 

The internal energy follows as U = F + TiSt + T 2 S 2 , 

U = ^JViTi + X -N 2 T 2 - gNTiT 2 (32) 

Notice that the sign of the order g correction is negative, 
due to entropic effects. To determine To from Eq. (||) we 
need entropies s[°^ , S 2 °^ , that can be read of from Eqs. 
©, @) at g = 0. One obtains T = T^T 2 2 and then 
from Eq. (§§) 

T g = T^T^[l + 2gv(2T^T^ -T x - T 2 )], (33) 

where v = N/(Ni + N 2 ), v k = N k /(N! + N 2 ), k = 1,2. 
According to Eq. (|l|) and using = N k T k /2, one gets 
for .g = the maximal amount of work 

W{0) = ~Nx (Ti - T? T2 2 ) + ^N 2 (T 2 - T» 1 T? ) (34) 
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Taking into account that V (T 1 ,T 2 ) = NTiT 2 , one gets 
from Eq. ( |23| ) a non-negative shift in the maximal of 
amount of work 



gWW=gN(T?T? 



Ti)(T 2 - T^T^ 2 ) > 0. (35) 



Only 



The equality is realized in the trivial case T% = T 2 . 
for g < this mechanism would reduce the work. 

Eq. (|2S|), and especially our model-dependent result 
(^5|), show that the full amount of the extracted work 
can increase due the direct coupling. Let us now return 
to the dissipative effects. For the considered model the 
energy dissipated per unit of time IT = T 2 S can be de- 
rived from Eq. (ph. It reads f| 



n 



8g^N_ 



7 r 2 (Ti - T 2 f + 0(j 2 ) 



(36) 



Our aim is now to obtain the characteristic time T, after 
which the energy dissipated due to the direct coupling is 
comparable with the energy gW^ gained according to 
Eq. (|35|). For the dissipated energy an upper estimate 
can be given as II T, and we get from II T = gW^\ 



8g 



T 2 (T^Tp - Ti)(T 2 - TpTp) 



(Ti - T 2 y 



(37) 



To be able to neglect the dissipated energy, the duration 
of our process t must be much smaller than T. On the 
other hand, since we are getting the work in the relax- 
ation process, its duration t must be much higher than 
the largest relaxation time T 2 . Thus, for times 



r 2 < t < r 



(38) 



it is possible to perform more work due to the presence 
of the direct coupling. The necessary condition T 2 -C T 
is realized mainly when g is small. If one was not able 
to complete the relaxation in the time- window (|38]), then 
for t ~ T the gained work will be equal to that without 
any direct coupling. The same analysis can be applied 
for the general case. 

So far we have compared the efficiencies of two systems 
with the same initial temperatures Ti, T 2 , and different 
values of g. One can also compare cases of identical initial 
energies, (Ti,T2,g_= 0) and (T\,T 2 ,g > 0), where the 
temperatures T\,T 2 are defined by U(Ti,T 2 ,g = Q) = 
U(Ti,T 2 , g > 0). The analysis is very similar to that 
given above, and indicates that our main result remains 
valid also in this case. There are examples of (Ti,T 2 ), 
for which the direct coupling enhances the work. 

In the context of our main result it is useful to inves- 
tigate which amount of work W(0 — > g) should be spent 
by external sources to switch on the small coupling g, 
starting from the state with g = 0. We shall consider the 
two extremes, very slow and very fast switching, which 
happen to give the same answer for small g. In the first 
case one uses Eqs. ( |l3| , [T^ ) to obtain (Tj , T 2 are constant) 



For the very fast switching the initial state does not 
change, and the main change comes from the Hamilto- 
nian (§): W/ Ost (0 g) = {H{g) - 7T(0)) ,which for 
small g leads to the same result as in Eq. (^9|). Using 
Eqs (|3[ |3^, [39]) one readily notices that there are tem- 
peratures, for which W a d(0 — > g) < gW^ l \ implying that 
the cost for the switching is less than the gain due to cou- 
pling: W(0) + W ad (0 ^g)< W(g). 

Conclusion. — Untill now it was believed that the pres- 
ence of a direct interaction between the baths of a heat 
engine reduces its efficiency The purpose of the 

present paper is to demonstrate that it can enhance the 
efficiency. Having changed the initial and final states, 
a direct coupling introduces, of course, both a change 
in work and dissipation. We show that the characteris- 
tic times of these two quantities can be well-separated. 
For times in the window (|3^) the work can be enhanced, 
though the dissipation is not yet relevant. This additional 
amount of work, which can be obtained from Eqs. ( |23| , 
p5| , [55j), is provided by the modified final state of the 
baths. 

Finally, we will briefly discuss related studies. 
Refs. |u]]ll|] consider the local thermodynamic efficiency 
of brownian motors and related models. The statement 
of this problem differs from the one considered by us, but 
it is interesting to mention that the role of a direct in- 
teraction between baths was studied also in this context 
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